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ABSTRACT 

Seven Class sources mapped with SCUBA at 850 and 450 /iin are modeled using a one dimensional 
radiative transfer code. The modeling takes into account heating from an internal protostar, heating 
from the ISRF, realistic beam effects, and chopping to model the normalized intensity profile and spectral 
energy distribution. Power law density models, n(r) oc r^^ , fit all of the sources; best fit values are mostly 
p = 1.8 ± 0.1, but two sources with aspherical emission contours have lower values {p ^ 1.1). Including 
all sources, {p) — 1.63 ± 0.33. Based on studies of the sensitivity of the best-fit p to variations in other 
input parameters, uncertainties in p for an envelope model are Ap — ±0.2. If an unresolved source (e.g., 
a disk) contributes 70% of the flux at the peak, p is lowered in this extreme case and Aj3 =tg j. The 
models allow a determination of the internal luminosity {{Lint) = 4.0 Lq) of the central protostar as well 
as a characteristic dust temperature for mass determination {{Tiso) — 13.8±2.4 K). We find that heating 
from the ISRF strongly affects the shape of the dust temperature profile and the normalized intensity 
profile, but does not contribute strongly to the overall bolometric luminosity of Class sources. There 
is little evidence for variation in the dust opacity as a function of distance from the central source. The 
data are well-fitted by dust opacities for coagulated dust grains with ice mantles (Ossenkopf & Henning 
1994). The density profile from an inside-out collapse model (Shu 1977) does not fit the data well, unless 
the infall radius is set so small as to make the density nearly a power-law. 

Subject headings: stars: formation — ISM: dust, extinction — ISM: clouds — ISM: individual (B335, 
B228, L723, IRAS03282-f 3035, L1527, L483, L1448C) 
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1. INTRODUCTION 

Modern theories of star formation predict the evolution 
of the density structure, n{f,t), and velocity structure, 
v{r,t), of the gas and dust envelope of protostars. Opti- 
cally thin dust emission at submillimeter wavelengths pro- 
vides an observational constraint on the density distribu- 
tion of the protostellar envelope and therefore constrains 
theoretical models of star formation. 

Class protostars represent an early, highly embedded 
phase during the formation of low mass stars (M < few 
Mq). The original evolutionary sequence for low mass 
protostars (Class I, II, and III) is based on the shape of 
the spectral energy distribution (SED) from observations 
at near and mid- infrared wavelengths (Lada 1987). Low 
mass protostars are thought to evolve from a thick dusty 
envelope where most of the energy is re-radiated by dust 
in the far-infrared (Class I) to progressively less embedded 
objects with near and mid-infrared excesses from dusty 
disks (Class II or classical T Tauri stars and Class III or 
weak-line T Tauri stars). The discovery of extremely em- 
bedded objects with submilfimeter telescopes led to a new 
class of protostars, Class objects, which are so highly 
enshrouded that their near infrared emission has generally 
not been detected and their SEDs peak longward of 100 



/im. Observationally, Class sources are cores that have 
Lboi/Lsmm < 200, where Lsmm is the total luminosity 
detected longward of 350 /im (Andre et al. 1993). Alter- 
natively, Class sources are characterized by Thoi < 70K, 
where Tboi is the temperature of a blackbody with the same 
mean frequency as the observed SED (Chen et al. 1995). 
Based on the relative numbers of Class and Class I ob- 
jects, it has been argued that the timescale for Class ob- 
jects is short, perhaps 10* years (Andre et al. 2000). Class 
objects have powerful outflows, which suggest high accre- 
tion rates (Bontemps et al. 1996). As data improve, some 
Class I sources are being reclassified to Class sources 
(e.g., Shirley et al. 2000, Young et al. 2001), suggesting 
a reexamination of this argument (see also Visser, Richer, 
& Chandler 2001). 

It is important to understand the structure of the en- 
velope of Class objects since they are likely to be the 
earliest observed phase of star formation with a central 
accreting protostar (Andre et al. 2000). The emission 
from an optically thin dust shell observed at an impact 
parameter, b, with dust opacity, k^, that does not vary 
with radius, is given by 



(Adams 1991) 
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where Tq is the outer radius. For an op- 
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tically thin envelope (at all wavelengths) dominated by a 
central source of luminosity, Lint, the dust temperature 
distribution can be approximated by a power law: 

Wcx(^)'^ a L«/,V«, (2) 

where q = 2/(4 + /?) (cf. Doty & Leung 1994) and /3 
is the power law exponent of the dust opacity (k oc z/''), 
which typically lies between 1 and 2 in the submillime- 
ter. If wc also assume the density distribution follows a 
power law, n(r) oc r~^, then in the Rayleigh- Jeans limit 
{Td » hv/k), the specific intensity integral simplifies to 
Iv oc 6^™, wlu^re rn = p + q — 1. Previous studios of the 
density structure of Class objects assumed a tempera- 
ture power law of the form Td{r) oc (Walker et al. 
1990, Ladd et al. 1991, Chandler & Richer 2000, Hoger- 
heijde & Sandell 2000, Shirley ct al. 2000, and Motte & 
Andre 2001). However, this approach is not valid in the 
outer envelopes of low luminosity sources due to a break- 
down in the Rayleigh- Jeans approximation at wavelengths 
shorter than 1 mm and due to external heating from the 
interstellar radiation field (ISRF) (Shirley et al. 2000). 
It then becomes necessary to calculate Td(r) self consis- 
tently in the integral in Equation (1) to reveal the density 
distribution of the envelope. 

In Paper I (Shirley ct al. 2000), 21 low mass cores within 
325 pc were observed at 850 iim and 450 i^m using the 
Submillimeter Common User Bolometer Array (SCUBA) 
(Holland et al. 1999) on the JCMT 15-m radio telescope. 
Thirteen sources from the original sample were classified 
as Class 0, three sources having been previously classified 
as Class I. In this paper, we present ID dust radiative 
transfer models for seven of the Class sources selected 
for spatial isolation and high signal-to-noise radial profiles 
i(^/^)peak = 60): B335, B228, L723, IRAS03282-F3035, 
L1448C, L1527, and L483. B335, B228, IRAS03282-F3035, 
and L1448C appear circular in the dust continuum maps 
down to the 20% contour, while L1527, L483, and L723 
are clearly aspherical at that level. The bolometric tem- 
perature ranges from 23K (IRAS03282+3035) to 52K 
(L483) and the bolometric luminosity ranges from 1.2 Lq 
(IRAS03282+3035 and B228) to 13 Lq (L483). In Paper 
II (Evans et al. 2001), three pre-protostellar cores (L1512, 
L1544, and L1689B) were modelled using one-dimensional 
radiative transfer and a beam convolution code. In this 
paper, the seven Class sources will be modeled using the 
same procedure used in Paper II, with the addition of an 
internal luminosity source. The modeling procedure and 
inputs are discussed in Section 2. We use B335 as a test 
object (Section 3.1) to model the sensitivity to the model 
input parameters. Individual sources are modeled in Sec- 
tions 3 and 4, while the implications of our best fit models 
are discussed in Section 5. 

2. Id dust radiative transfer modeling 

The ID dust models are calculated from a modified ver- 
sion of the Egan, Leung, and Spagna (1984) continuum ra- 
diative transfer code. This code iteratively calculates the 
equilibrium dust temperature on a ID radial grid by simul- 
taneously solving the combined moment radiative trans- 
port equations in quasi- diffusion form and the energy bal- 
ance equations as a two-point boundary value problem. 



The radiation field is constructed by solving a sot of ray 
equations along impact parameters, b, through the cloud. 
There are four physical inputs to the radiative transfer 
code: the density distribution, n(r); the internal source 
luminosity. Lint] the scaling factor for the interstellar ra- 
diation field, sisRF', and the dust opacity, k,j. In addi- 
tion, the radial grid (100 points) and impact parameters, 
as well as the wavelength grid (59 wavelengths) are chosen 
to cover the relevant range so that results are insensitive 
to the details of these grids. The equilibrium dust tem- 
perature distribution, Id(r), which is the output from the 
dust code, is used by an observation simulation code (Pa- 
per II) to calculate the normalized radial intensity profiles, 
I'^°™{b), by solving Equation (1), performing beam con- 
volutions, and simulating chopping. The code also models 
the observed SED, Sv{Omb), by convolving the model in- 
tensity with the beam {6mb) used in each observation. 

The internal radiation field is assumed to be a black- 
body with an effective temperature of 6000 K. Since all 
of the objects we are modeling are very opaque at near- 
infrared wavelengths, the exact shape of the spectrum of 
the internal radiation field is not important; only the total 
internal luminosity is important. The shape of the ISRF 
was determined from COBE results (Black 1994) plus the 
cosmic microwave background (CMB) and the ultraviolet 
component of the ISRF (A < 0.36/xm) from van Dishoeck 
(1988). This ISRF is considerably stronger in the infrared 
than the previous versions of the ISRF (Mathis, Mezgcr, 
& Panagia 1983) (see Figure 2 of Paper II for a plot of 
the different ISRFs). We modify the strength of the ISRF 
by multiplying all portions of the ISRF spectrum except 
the CMB with a factor denoted by sisrf- For the mod- 
els discussed in this paper, we used a coarse grid with 
Sisrf = 0.3, 1.0, and 3.0. These factors for the far-UV 
(FUV) portion of the ISRF correspond to 0.45Go, 1.5Go, 
and 4.5Go respectively, where Go is integrated FUV flux 
between 91.2 nm and 220 nm in units of 1.6 x 10^ erg s~^ 
cm"2 (^(,f Hollenbach et al. 1991). In Paper II, the best 
fit to the observed radial profiles and SED suggest a lower 
strength to the ISRF {sisrf ~ 0.3 - 0.5). 

The dust opacity, , was taken from the models of Os- 
senkopf and Henning (1994) for grains that have coagu- 
lated for 10'' yr at a density of 10^ cm~^, both with (0H5) 
and without (0H2) accreted ice mantles. We assume a 
gas to dust ratio of 100 by mass to convert 0H5 opaci- 
ties (per gram of dust) to opacity per gram of gas. 0H5 
opacities have been successful at reproducing the SED of 
both high mass (van der Tak et al. 2000) and low mass 
star forming cores (Paper II). 0H5 opacities at submil- 
limeter wavelengths can be approximated by a power law 
[ku oc u^) with (} ~ 1.85 and an opacity of 1.8 x 10^^ 
cm^ per gram of gas at 850 jiui. The 0H2 dust opacities 
are lower at far-infrared wavelengths but up to 1.6 times 
higher at submillimeter and millimeter wavelengths than 
0H5 dust opacities. The cross-over point is near 350 /xm 
(see Fig. 2 of Paper II). 

The JCMT beams derived from AFGL618 and Uranus 
were used for the beam convolution of the model profiles 
at 850 and 450 /im (see Figure 1). It is very important 
to convolve the model specific intensity distribution with 
a realistic beam profile to derive an accurate estimate of 
the envelope density distribution. Using a gaussian beam 
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shape instead of the actual beam profile can profoundly 
distort the interpretation of the power law index, p. by up 
to 0.5 (Shirley et al. 2000). Since the beam shape was 
stable during the second half of the night in April 1998 
(see Paper I), nine Uranus profiles were added together to 
use an average beam profile with high signal-to-noise out 
to 70" from the center of the map. No planets were avail- 
able to make beam maps for the January 1998 sources 
(L1527, IRAS03282, L1448C); so the secondary calibra- 
tor, AFGL618 was observed. The January 1998 observa- 
tions were taken in the first half of the night when the 
beam shape is not so stable and changes shape continu- 
ously as the telescope cools. AFGL618 is much weaker 
than Uranus at 850 and 450 /im; therefore the beam maps 
have much lower signal-to-noise at large radii. The Jan- 
uary beam map was produced by averaging together three 
AFGL618 maps observed during the same time of night 
that our objects were observed. Because the 450 ^m pro- 
file of AFGL618 becomes too noisy beyond 34", the aver- 
age Uranus beam profile from April was used for the 450 
/xm January beam beyond 34". These beam calibration 
difficulties introduce larger uncertainties, ~ 0.1, in the 
best fit models for the January sources. 

Many of the radial profiles from Paper I show a turn- 
down in the intensity profile beyond 60" from the center 
of the map. This turn-down could be due to a steepening 
of the density profile or due to the effects of chopping. To 
test the possibility of a steeper density profile, we must ac- 
count for the chop throw. The SCUBA observations were 
chopped in azimuth with a 120" chop throw. Since the 
SCUBA array is 2D, all of the positions within a single an- 
nulus chop different distances from the center of the map. 
Therefore, our ID model can only include an approximate 
simulation of the true effects of chopping. The details of 
how we simulate chopping has a noticeable effect on the 
shape of the model intensity profiles beyond 60"; conse- 
quently, we do not attempt to model normalized intensity 
profiles beyond this radius. 

The agreement between the model and the data can be 
quantified in terms of the reduced chi-squared for the nor- 
malized intensity profile at wavelength A: 

■ / rnormfi. \\ { jnormfi. \\ ") 2 
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I Nb , (3) 



where I^°™{bi) is the azimuthally averaged, normalized 
intensity in a circular aperture at impact parameter 6,, 
(Jiihi) is the uncertainty in the data, and is the num- 
ber of impact parameters. Only points spaced by a full 
beam are used in computing x\ to avoid introducing cor- 
relations. We calculate x| for the 850 and 450 /im profiles 
and define as the sum of x^r^Q and xlso- The signal-to- 
noise was higher for the 850 /um maps and therefore the 
xiso has the most weight in determining the best fit. 

The agreement between the model SED and the ob- 
served one is quantified by Xsed ' calculated from a similar 
equation: 

cobs Qmod ^ 



Xsed — 



(4) 



where 5'°^'* is the observed flux into a beam and 5'™°'' is 
the modeled flux into the same beam. When photometry 
at the same wavelength with different beams exists, the 



points are both considered in calculating x%ed- Wave- 
lengths shorter than 60 /im are not included bec;ause they 
are expected to be optically thick and therefore very sensi- 
tive to asymmetric geometries (e.g., outflow cavities, flat- 
tened envelopes), which we are currently unable to model. 
In addition, much of the far-infrared photometry has large 
and uncertain calibration errorbars, and the opacities as a 
function of frequency are uncertain. For all these reasons, 
poor fits to the SED are not considered a serious problem 
in constraining the density distribution; however we often 
find that the best Xri which considers only the profiles, 
occurs for the model with the best xLt = Xr + Xsed- 

In fact, the SED and the normalized radial intensity 
profiles provide nearly orthogonal constraints on model 
parameters (§3). The SED is sensitive to the strength of 
the ISRF and a mass x opacity product (Paper II). In par- 
ticular, the flux density at 850 /im constrains the mass, 
while the full SED provides information on the the varia- 
tion of dust opacity with frequency, subject to the caveats 
mentioned above. Measurements of the flux density into 
different beams at the same wavelength provides some con- 
straint on the density distribution, but the shape of the 
density distribution is much better constrained by the nor- 
malized radial profile {I"°^"^{b)). 

The simplest models to test are single power laws. 



n(r) 



r G [ri,ro] 



(5) 



in cm~'^ for the gas density. The density, n/, is normalized 
to a fiducial radius, r/, of 1000 AU and represents the total 
gas number density (n = n(H2) + n(He) -|- ... = p/ (/xmn), 
/i = 2.29). There are seven parameters in the power law 
models: the power law exponent, p: the density at a fidu- 
cial radius, n/; the inner radius, r^; the outer radius, rg', 
the internal source luminosity. Lint; the dust opacity, 
and the strength of the interstellar radiation field, sjsrf- 
In these models, the shape of the density profile, defined 
by p is constrained by I"°™{b), and only slightly affected 
by other parameters. In contrast, n/ and hence the mass, 
are constrained by the observed flux density at optically 
thin wavelengths; we use 5*^(850) because the calibration 
errors are lowest. The resulting mass depends inversely 
on the opacity at 850 /xm {k„{850)) and weakly on other 
parameters. 

We also test models of inside out collapse (Shu 1977), 

hereafter referred to as Shu77 models. These are char- 
acterized by seven parameters: Ce//, the effective sound 
speed; the radius that encloses the infalling gas, ri„/; rf, 
fo] Lint; K„; and sjsrf- Inside Tin/, the density distribu- 
tion tends toward n(r) oc r~^'^, while n{r) oc outside 
Vinf ■ The shape of the density profile is thus set by Vinf , 
while the normalization is set by fle//- We fix fle// based 
on observations of optically thin spectral lines; thus Shu77 
models have no freedom in the normalization of the density 
profile, unlike power law models. 

The constraints on the other parameters arc similar for 
the two types of density profiles. The internal source lumi- 
nosity (Lint) is constrained by the integral of the SED and 
secondarily by sisrf- The internal luminosity dominates 
the bolometric luminosity over heating from the ISRF for 
sources with Lint ^ IL©. The model internal luminos- 
ity is tuned until the model bolometric luminosity {L^'^'^) 
matches the observed bolometric luminosity {Lobs, Paper 
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I), using the same method to integrate over the SED. The 
model bolometric himinosity is c;alc;ulated using the same 
wavelengths and beams as the observed bolometric lumi- 
nosity. Lobs and L'^"'' may be greater than Lint for low 
luminosity sources because the ISRF adds energy. In some 
cases, especially for more luminous sources Lobs underes- 
timates Lini because the beams used at some wavelengths 
do not capture all the emission (see Butner et al. 1990 for 
a full discussion). The choice of the overall opacity law 
{ki,{i')) is constrained by the shape of the SED once other 
parameters, like p, Lint, and sisrf are constrained. 

Models of B335 are used to test our assumptions of 
nearly orthogonal constraints and the effects of changing 
various parameters in the model (§3). 

3. TESTING MODEL PARAMETERS B335 

B335 (IRAS 19345-F0727, L663) is an extensively stud- 
ied Class object within the Barnard 335 dark cloud. Its 
submillimeter emission is very nearly circularly symmet- 
ric (Huard et al. 1999, Paper I, Motte et al. 2001). The 
core is the one of the best cases for a collapse candidate as 
deduced from models of observed asymmetric line profiles 
in CS and H2CO (Zhou et al. 1993, Choi et al. 1995). 
Rotation, if any, is very slow (Frerking et al. 1987, Zhou, 
1995), making simple spherical models reasonable. B335 
has an outflow that lies nearly in the plane of the sky along 
an east-west direction (Goldsmith et al. 1984). There is 
little direct evidence in the submillimeter continuum maps 
of extensions or flattening along or perpendicular to the 
outflow direction, making this core a suitable choice for 
ID modeling. Harvey et al. (2001) have recently studied 
B335 using a near-infrared extinction mapping technique 
to probe the density structure (cf. Alves, Lada, & Lada 
2001). These observations provide an important check on 
the consistency of our submillimeter continuum models. 
We shall model the density structure of the outer envelope 
of B335 using both a single power law and the inside-out 
collapse model of Shu (1977). 

3.1. Power Law Models 

A p = 2 power law is a good starting model because it 
has the same density distribution as a singular isothermal 
sphere, until it is truncated at an inner radius. The outer 
radius was chosen to allow proper simulation of chopping 
at the distance of B335 (250 pc; Tomita et al. 1979). The 
chop throw for our observations was 120", corresponding 
to 30,000 AU at the distance of B335. Both the radial 
proflle and the observed beam profile contain useful in- 
formation 60" from the center of the map; therefore, to 
allow the beam to chop onto the model density distribu- 
tion for radial points 60" from the center, the modelled 
density distribution must extend to twice the chop throw 
(240") or 60,000 AU at the distance of B335. NICMOS ob- 
servations of B335 indicate that reddening of background 
stars blend into the background noise beyond 125" (Har- 
vey et al. 2001), which is within the range of the outer 
radii tested in our models. We initially chose the inner 
radius such that ro/rj = 1000 so that rj = 60AU for B335. 
For the initial test model, we used 0H5 opacities. The 
fiducial density, n/, was varied to match the observed flux 
at 850 microns; S'gfg = 3.91 ± 0.22, and the internal lu- 
minosity, Lint, was varied until the model Lhoi matched 



Lobs = 3.1 ± O.IL0 (Paper I). The resulting best fit was 
Hf = 1.7 X lO^cm^"* and Lint = 2.5L0 (see Table 1). This 
basic model was used to test the eflFects of varying Vo, Vi, 
Ki/, p, Sisrf, and distance. The sensitivity of p and Lint 
to changes in the other parameters is summarized in Table 
2. 

Changes in ro result in changes in the total envelope 
mass proportional to r^~P for a power law; however, the 
column density along a line of sight, and hence the emis- 
sion at long wavelengths into a fixed beam becomes insen- 
sitive to To for p > 1 {N oc rl~^{l — {ri/roY~^)- For the 
test model, the density is 4.7 x 10^ cm'^ at 60,000 AU. 
This low value is probably not realistic as the dense core 
is embedded in an extended cloud (Frerking et al. 1987). 
Even if a constant density of 10^ cm^'^ is used for the lower 
limit to the density in the envelope, doubling the outer ra- 
dius to 120,000 AU only increases Ay by slightly less than 
1 magnitude of visual extinction. This lower limit is used 
for all subsequent models. Allowing the density to fall be- 
low 10'^ cm~^ results in no significant change in the best 
fit fiducial density or internal luminosity. Factor of two 
changes in To have negligible effect on the best fit param- 
eters. If To is made small enough {ro < 15000 AU), the fit 
degrades (Table 1), but this is smaller than the extent 
of observed emission in the SCUBA map. 

While the column density increases as rj decreases, 
beam dilution is more important. The SCUBA beams 
are much larger than the inner radii in our models (the 
450 /zm beamwidth of 8" corresponds to a radius of 520 
AU for our nearest source, B228); therefore changes in the 
inner radius will not have substantial effects unless the in- 
ner radius gets too large. Factor of two changes in have 
negligible effects; a factor of 4 increase in begins to de- 
grade the fit to the 60 /um point on the SED (Figure 2). 
However, 60 /zm observations can be affected by deviations 
from spherical geometry, so the 60 /im fiux should not be 
considered a strong constraint on one dimensional models. 
We do not model the effects of a disk (see Section 5.5). 

The result of these tests is that the best-fit model pa- 
rameters are not sensitive to changes of factors of 2 in Vo 
or ri. Larger changes worsen the fit, indicating that there 
is no evidence in our data for either inner or outer bound- 
aries in the power laws. 

Calculations of dust opacities differ substantially be- 
tween dust models, especially at long wavelengths; for 
example 0H5 and 0H2 dust models have opacities at 
A = 1000 /um that are 6 and 12 times those of the model 
for dust in the diffuse interstellar medium (Draine & Lee 
1984). With no further information, these differences 
would result in substantial uncertainties in rif, or equiv- 
alently the mass since S^, cx Mk,^. However, the SED is 
strongly affected by the choice of opacity law and can be 
used to distinguish dust models (e.g., Butner et al. 1991, 
van dcr Tak et al. 2000) if the density distribution is inde- 
pendently constrained. If 0H2 opacities are used instead 
of 0H5 opacities, the best fit power law remains p = 2.0 
with n f about one-third that of the test model. The higher 
K,y(850) of 0H2 dust requires less column density to match 
the flux at 850 /im. None of the power law flts using 0H2 
dust produce enough emission to match the SED at 100 
/zm > A > 350 fxm (Figure 2), yielding quite large values of 
XsED- Changing the submillimeter opacity by a factor of a 



5 



few does not strongly affect the shape of the model density 
profiles, confirming our conclusions about the orthogonal- 
ity of the radial profile and the SED. A dust opacity that 
varies with radius could in principle affect the fit to the 
shape of the density profile: however, we see no evidence 
for variation of opacity with radius (see §5.3). 

Other models of star formation predict density distri- 
butions different from the singular isothermal sphere (e.g. 
McLaughlin & Pudritz 1997, Ciolek & Basu 2000). With- 
out explaining the details of these models we note that the 
differences often can be approximated by variations in the 
power law index (p) in the bulk of the envelope. We there- 
fore considered models with p = 1.5, 2.0, and 2.5 (Figure 
3). Both the p = 1.5 and p = 2.5 models fit very poorly 
(Fig. 3). The p = 2.5 model shows that the predicted ra- 
dial profile becomes very sensitive to the beam shape for 
steep power laws. Variations in p of 0.1 produce noticeable 
changes in the intensity profiles. In summary, we find that 
the p = 2.0 initial test model is the best fit profile for 0H5 
dust and sjsrf = 1-0, and that variations in p of 0.1 can 
be distinguished if all other parameters are fixed. 

Heating by the ISRF affects the temperature distribu- 
tion in the outer envelope. The temperature levels off and 
eventually rises again towards the outer boundary of the 
cloud. The location of the minimum in the temperature is 
dependent on the ISRF strength and Lj„t (Figure 4). The 
strength of the ISRF could affect the intensity profiles and 
the SED. We tested these effects on the best fit power law 
by multiplying the ultraviolet to far-infrared portion of 
the ISRF by factors of sjsrf = 3 and sjsrf = 0.3. With 
SiSRF = 3.0, the best fit is obtained with p = 2.2 and 
Lint = 1-5 Lq. Obviously, the greater contribution to the 
flux from the ISRF explains why a lower internal lumi- 
nosity is appropriate. Even though Xr and Xsed are not 
significantly different from the standard model, the pre- 
dicted SED is too bright at FIR wavelengths and severely 
underestimates the 60 /im flux. Using sjsrf = 0.3, the 
best fit is obtained with p = 1.8 and Lint = 3.3L(5 (Fig- 
ure 4 & 5). This model clearly yields the best fit to the 
intensity profiles of all of the power law models that we 
have considered, but the model SED has too much flux 
at shorter wavelengths (i.e., 60 /um). A lower strength of 
the ISRF is consistent with the best fit models in Paper 
II. A study of [CII] and [01] emission lines with ISO to- 
wards B335 is consistent with a slightly higher strength 
of the ISRF, 2 - 3Go (Nisini et al. 1999). Clearly, there 
is uncertainty in the strength of the ISRF around B335. 
We conclude that uncertainties of a factor of 10 in sjsrf 
cause changes in p by ±0.2, but that the SED provides 
some constraints that can decrease this uncertainty. 

The effect of the uncertainty in distance was checked on 
the initial test power law by assuming B335 was at half 
the published distance (125 pc; cf. Harvey et al. 2001). 
Since the core is now closer, smaller inner and outer radii 
(30 pc and 30,000 pc) were used. The best fit power law 
index and shape of the model profile are very insensitive to 
the uncertainty in the distance {Xr increases by less than 
3%). The best fit ij„t = 0.7 Lq, and the total mass in the 
model decreased by a factor of 2. The mass of a power law 
envelope scales as Menv oc Ufr^^P. For factor of 2 uncer- 
tainties in the distance, Lint will be uncertain by factors 
of about 4, the mass by factors of 2 (for p = 2), but the 



best fit p is unaffected. 

The effect of changing the observed beam shape was 
tested by substituting the January 1998 beamshape into 
the best fit power law model for B335. The difference in 
beamshape between January and April 1998 beams was 
greater than the difference between individual beam pro- 
files observed during April 1998. The best fit power law 
was slightly lower, p = 1.7, indicating that uncertainties in 
the beam shape may effect the interpretation of the best 
fit model by as much as Ap ^0.1. 

A p — 1.8 power law with sisrf = 0.3 provides the best 
fit to the radial profiles for B335 (Figure 5). However, the 
Xsed ^ good for this model as for the test model. 

A somewhat different opacity law would fit better. The fit 
to the radial profiles is not perfect, with some systematic 
deviations as a function of impact parameter that suggest 
deviations from a simple power law, but given uncertain- 
ties in the beam shape, we consider this fit to be adequate. 

In summary, our tests support the idea that the normal- 
ized radial profile (/"°''™(&)) and the SED provide largely 
orthogonal constraints. 7"°'''" (6) constrains the shape {p) 
of the density profile, while 51,(850) constrains the nor- 
malizing density, n/, and hence the mass. The density 
and mass enter as a product with the opacity at 850 fim, 
leading to uncertainties in the mass equal to the uncer- 
tainties in the opacities. However, the rest of the SED 
depends on the overall opacity law; within the limitations 
of spherical models, the overall SED constrains the opacity 
law and hence reduces the uncertainty in mass. Variations 
in other parameters, such as opacity, distance, and Lint, 
have negligible effect on p (Table 2). Variations in the 
ISRF have the largest effects: a factor of 3 in either direc- 
tion in Sisrf cause variations in p of ±0.2. Factors of 3 in 
Sisrf can be constrained by the SED, so larger errors are 
unlikely. Opacity and the strength of the ISRF strongly 
affect Lint (Table 2). 

3.2. Shu77 Inside- Out Collapse Models 

The Shu77 density distribution has successfully matched 
asymmetric line profiles seen towards Class sources in- 
cluding B335 (Zhou et al. 1993, Choi et al. 1995, Hoger- 
heijde et al. 2000). Based on Monte Carlo radiative trans- 
fer models of CS and H2CO asymmetric line profiles, Choi 
et al. found that a Shu77 model with r,„/ = 6200 AU 
and Ueff = 0.23 km s~^ was the best fit for an inside-out 
collapse model for B335. We mostly constrain a^ff to this 
value, unlike the procedure used by Hogerheijde & Sandell 
(2000), who allowed Oe// to be a free parameter. This lack 
of flexibility in the normalizing density puts Shu77 models 
at a disadvantage compared to the power law models with 
two free parameters. 

The Shu77 density distribution with rinfaii = 6200 AU 
was tested assuming sjsrf — 1-0, and 0H5 opacities. 
Compared to the best-fit power law model, the density 
in this Shu77 model is about 5 times less. There is sim- 
ply not enough material to match the SED: Lint had to 
be increased to 6.5 Lq to match Lobs because the lower 
optical depths allowed too much radiation at A < 60 /^m, 
which is not observed. Even with the large Lint, the fits 
to both the SED and radial profiles are extremely poor 
(Figure 6). The model radial profiles are too flat. This is 
not surprising since a p = 1.5 power law model did not fit 
the observed radial profiles. 
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Shu77 models would need smaller infall radii to pre- 
serve an r^^ density distribution throughout a greater ex- 
tent of the envelope. Smaller infall radii provide better 
fits to the radial profile. The best of these models, with 
Tinf = 1000 AU and sisrf = 0.3, has a very good Xr- 
However, r,„/ = 1000 AU is less than the FWHM of the 
SCUBA beam at 450 y^m. Thus, the best fit Shu77 model 
actually resembles a p = 2 power law in all portions of the 
profile except the central beam. A small infall radius is 
highly unfavored by models of profiles of molecular lines 
(x^ is 20 times worse, Choi ct al 1995). Further radiative 
transfer modeling of the molecular line emission using the 
best fit density profiles from the dust models and consider- 
ing possible abundance variations in the molecular tracers 
may resolve this discrepancy. 

Thus, the Choi et al. model has two problems: the den- 
sity is too low to match the observed SED, and the density 
distribution is not steep enough to match I™™{b). Be- 
cause other opacity models combined with Shu77 models 
have fit the SED of B335 (Zhou et al. 1990), the first 
problem is not necessarily fatal. The bigger problem is 
the failure to match I^'"-"'{b). 

The near infrared extinction study of Harvey et al. 
(2001) found that the Choi et al. (1995) parameters for a 
Shu77 model did fit the density structure of the outer en- 
velope, but require a scaling of the density by a factor of 5. 
Such a large increase would require a higher Oe/ / or B335 
to be much closer. A higher a^ff is not supported by the 
narrow linewidth observed towards B335 (e.g., Mardones 
et al. 1997). As shown in the power law models of B335, a 
closer distance does not strongly affect the interpretation 
of the best fit radial profiles; therefore, an infall radius of 
6200 AU would still not fit the observed normalized in- 
tensity profile. This solution, which worked for Harvey 
et al., does not work for our data. A small infall radius 
would be required to match the SCUBA profiles with a 
Shu77 model. The increased density required by Harvey 
ct al. (2001) to model their near-infrared extinction mea- 
surements makes their densities agree well with our best- 
fit power law model, which is based on 0H5 opacities at 
850 /im. The 0H5 opacity at 850 /im and the extinction 
used by Harvey et al. in the near-infrared thus produce 
consistent density estimates. This agreement supports the 
validity of 0H5 dust, as does our comparison of dust and 
virial masses (§5.2). 

4. SOURCES 

The other sources were modeled with less extensive ex- 
ploration of parameter space, guided by the results of the 
B335 models. 

4.1. B228 

B228 is a deeply embedded protostar (IRAS 15398-3359) 
which appears very similar to B335 at 850 and 450 /xm. 
The B228 contours are nearly circularly symmetric with 
extended envelope emission and a high signal-to-noise pro- 
file. Unlike B335, B228 is much less studied and therefore 
has less published SED information (see Paper I). B228 is 
an excellent candidate for ID dust modeling of the enve- 
lope emission and should be the subject of further detailed 
study. We shall try both power law and Shu77 models. 

A range of power law indices were tried with noticeable 



changes in the Xr with changes of Ap = 0.1. The best fit 
power law is p = 1.9 with a lower ISRF {sisrf = 0.3) and 
an internal luminosity of 1.0 Lq (see Figure 7). Since B228 
is closer than B335 (130 pc, see Paper I), the outer radius 
was decreased to 30,000 AU. While the p = 1.9 power law 
is a reasonable fit to both the profile and SED, it has the 
same problem as the B335 power law fits. The model pro- 
file is too steep in the inner portion of the envelope and too 
shallow in the outer portion of the profile. The deviations 
may be due to subtle changes in the beam shape or to de- 
partures from a single power law. Far-infrared photometry 
with better spatial resolution and more wavelength points 
would help to constrain the SED. 

The best fit Shu77 collapse model has an infall radius of 
1000 AU and an effective sound speed of 0.23 km s~^. The 
linewidth of optically thin lines in B228 (e.g., Au(N2H+) 
in Mardones et al 1997) is similar to that of lines in B335. 
Therefore an Oe// = 0.23 km s~^ is reasonable. As with 
B335, the infall radius is within the central SCUBA beam. 
The best fit Shu77 model was able to simultaneously match 
the observed bolometric luminosity and the fiux at 850 /zm. 
The fit to the profiles are not as good as the fit of the best 
fit power law, and Xsed considerably worse. 

4.2. L723 

The submillimeter emission from L723 (IRAS 
19156-1-1906) has asymmetries only in the lowest contours 
of the 850 and 450 fim emission. The higher contours 
are very circularly symmetric. A quadrupolar outfiow 
has been detected, indicating that L723 may contain a 
close binary (Palacios & Eiroa, 1999) within a common 
envelope. Alternatively, the outflow may be due to a sin- 
gle source with a very large opening angle (Hirano et al. 
1998). The largest extensions in the lowest contours of 
the SCUBA maps are roughly aligned with an east-west 
outflow. The unique characteristics of the outflow struc- 
ture and the possibility that L723 may be a proto-binary 
system warrant a detailed study of the dust continuum 
structure of the envelope. 

L723 is best fitted by a p = 1.8 power law with sjsrf = 
0.3 (Figure 8). The fit to the profiles and SED is very good 
using 0H5 dust. The model internal luminosity is 2.6 Lq 
indicating a contribution of 0.5 Lq to the bolometric lu- 
minosity from the ISRF. The large contribution from the 
ISRF may be due to the large size of model (60, 000 AU) 
and the greater distance (300 pc; Goldsmith et al. 1984) of 
L723. Emission was detected beyond 65" (-- 20,000 AU) 
in the 850/Ltm SCUBA map, larger than the outer radius 
quoted by Motte & Andre (2001) of 14,000 AU from 1.3 
mm continuum maps. As seen in the test models of B335, 
decreasing the outer radius to 30, 000 AU has a negligible 
effect on the shape of the radial profile and the interpre- 
tation of the best fit power law. 

The best fit power law is very similar to the other cir- 
cularly symmetric cores B335 and B228. While L723 may 
be a proto-binary system, there is no strong evidence for 
differences in the one dimensional dust models of L723 
and B335 or B228. Multiple dimensional dust models are 
needed to probe the extended structure at low contour lev- 
els in the SCUBA maps. The envelope structure of L723 is 
complicated by the extensive observed outflows. The im- 
pact the outflow structure has on the overall distribution 
of material in the envelope cannot be effectively modeled 
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with only one spatial dimension. 

The best fit Slni77 model has a small infall radius (1000 
AU) and effective sound speed of 0.29 km s~^. The fit to 
the profile and SED is very good. As with B335 and B228, 
the infall radius is within the central 450 /im SCUBA 
beam. There is no evidence for a break in the power law 
density distribution in the outer envelope beyond the cen- 
tral beam. 

4.3. IRAS03282+3035 

IRAS 03282+3035 is a nearly circular symmetric core 
with slight asymmetries in the lowest contour of the 850 
and 450 /um emission. There is a well studied outflow 
(e.g., Bachiller et al. 1994) that is nearly perpendicular 
to the extended submillimeter continuum emission. The 
signal-to-noise of the radial profiles was lower for this ob- 
ject (Paper I). 

The best fit to IRAS 03282+3035 is a p = 1.9 power law 
with lower strength of the ISRF, sjsrf = 0.3 (see Figure 
9). The model fits the 850 /im profile marginally better 
than the 450 fim profile. Since IRAS 03282+3035' was ob- 
served during January 1998, the beam profile is much more 
uncertain, resulting in larger uncertainties in the best fit 
model. Lint = 1-0 Lq for the best fit model, accounting 
for nearly all of the bolometric luminosity. 

We continue to use a distance of 220 pc (Ccrnis 1990) 
for consistency with Paper I, but the recent study of de 
Zeeuw et al. (1999) finds a greater distance of 318 pc us- 
ing Hipparcos parallaxes of the Perseus OB association. 
At the larger distance, the internal luminosity and mass 
in the envelope would increase by a factor of two. As was 
seen in the models of B335, the best fit power law index 
is not sensitive to distance and remains p= 1.9. 

The best fit Shu77 model again has a infall radius of 1000 
AU within the central SCUBA beam. Molecular line ob- 
servations towards IRAS03282+3035 indicate linewidths 
that are similar to B228 and B335 (A?;(N2H+) = 0.49 
km s~^, Mardones et al. 1997); however, a higher effective 
sound speed (0.26 km s"-'^) than used for B228 is needed 
to match the 850 /xm fiux and Lboi simultaneously using 
0H5 dust. 

4.4. L1U8C 

L1448C (also called L1448-mm) is a well studied Class 

protostar in the vicinity of 3 other deeply embedded 
protostars in the Perseus molecular cloud (L1448N(A & 
B), L1448NW; see O'Linger et al. 1999, Barsony et al. 
1998). L1448C drives a powerful, highly coUimated out- 
flow (Bachiller et al. 1990) for which proper motion has 
been detected in extremely high velocity SiO maps (Girat 
& Acord 2001). H2O mascr emission has been detected 
(Chernin 1995) indicating the presence of very dense gas. 
Recently, the infrared spectrum from 6/im to 190/xm has 
been observed with ISO and extensively studied (Giannini 
et al. 2001, Nisini et al. 1999). As for IRAS03282+3035, 
we continue to use a distance of 220 pc, but 318 pc is more 
likely (de Zeeuw et al. 1999); the effects on the model pa- 
rameters are described in the previous section. 

The submillimeter contours are circularly symmetric 
with a weak bridge of emission extending between L1448C 
and L1448N/L1448NW (Barsony et al. 1998, Chandler & 
Richer 2000, Paper I). L1448C is located 82" from the 



nearest source (L1448N); therefore we shall attempt to 
model L1448C as an isolated source with the strong caveat 
that the outer portion of the observed radial profile (i.e., 
> 45") is contaminated with emission from L1448N and 
L1448NW. Only points in the radial profile < 45" from 
L1448C are modeled. 

The best fit power law was p = 1.7 with sjsrf = 1-0 
(Figure 10). An outer radius of 45,000 AU was used for 
the best fit model. This model cannot be compared di- 
rectly to observations beyond half the distance between 
L1448C and L1448N/L1448NW (9000 AU projected on 
the sky). However, the large outer radius is appropriate 
for the model since all of the sources in the SCUBA map 
are clearly embedded in a diffuse envelope that extends to 
the edge of the map. The model of the 450 profile falls 
off too steeply for the best-fit 850 /um profile, which could 
be caused by a larger sisrf or uncertainties in the beam. 

This power law density structure is slightly flatter than 
that of the other protostars with symmetric submillimeter 
emission. If the entire profile of L1448C is contaminated 
from emission from L1448N/L1448NW, then the observed 
intensity profile would be flatter than for an isolated source 
resulting in a flatter best fit p. However, it seems unlikely 
that the northern sources could significantly alter the con- 
clusions of the modeling of the entire profile since they are 
located at least 18000 AU away. 

This is the only modeled source for which sjsrf = 1 
gave the best fit. This source is also the only source mod- 
eled with multiple cores observed within the SCUBA map. 
It is possible that this indicates that the ISRF is indeed 
greater in the vicinity of L1448C; however, the overall 
change in the Xtot between sjsrf = 0.3 and sjsrf = 1.0 
is only 25%. 

4.5. L1527 

L1527 is a well studied Class source (IRAS 
04368+2557) in the Taurus molecular cloud complex char- 
acterized by asymmetric submillimeter and far-infrared 
emission extending from the southeast to northwest direc- 
tion (Ladd et al. 1991, Chandler & Richer 2000, Hoger- 
heijde & Sandell 2000, Paper I, Motte & Andre 2001), 
strong evidence for rotation (Goodman et al. 1993, Zhou 
et al. 1996, Ohashi et al 1997), a molecular outflow in the 
east-west direction (MacLeod et al. 1994, Bontemps et 
al. 1996), and an associated near-infrared nebula (Eiroa 
et al. 1994). The normalized radial intensity profile was 
not well fitted by a single power law in Paper I, as it shows 
considerable curvature. Previous studies of the dust con- 
tinuum emission (Hogerheijde & Sandell 2000, Motte & 
Andre 2001) have estimated the density power law to be 
near p = 1, much lower than the best fit power laws of 
B335 and B228. 

The best fit power law model (see Figure 11) is a shal- 
low p = 1.1 power law with sjsrf = 0.3, in agreement 
with previous dust continuum studies of L1527. The fit to 
the 850 /im profile is within the errorbars except for the 
last point. The observed profile is slightly fiatter than the 
model p = 1.1 power law between 5000 and 6500 AU. The 
450 /zm model profile is fiatter than the observed profile in 
the inner portion but matches extremely well in the outer 
envelope. These discrepancies may be explained by the in- 
ability of the one-dimensional dust code to appropriately 
model the observed large asymmetries in the dust contin- 
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uum maps or may be explained by the problematic Jan- 
uary 1998 beam profiles (cf. Section 2). The beam effects 
cannot be strong enough to change our conclusion that the 
density structure in LI 527 is clearly flatter than the den- 
sity structure of B335 or B228. The curvature in I^°™{h) 
can be matched by simple power laws when Td{r) and the 
effects of the beam and chopping are properly sinndated. 
The fit to the SED is poor at far-infrared wavelength as 
indicated by the the poor xIed for all of the models in 
Table 1. 

The sensitivity of the model to asymmetries can be 
tested by eliminating sectors in the azimuthally averaged 
intensity profile. L1527 is elongated in a south-east to 
north-west direction. If a symmetrical sector of 70° cen- 
tered along the major axis is removed from the azimuthal 
average, then the best fit power law index changes only to 
p = 1.0. Therefore, the density distribution appears to fall 
off more slowly even perpendicular to the extension; thus 
the lower value of p appears to be real and not an arti- 
fact of azimuthally averaging an elongated intensity map. 
Uncertainties on the order of Ap ~ ±0.1 result from az- 
imuthally averaging the L1527 continuum emission. The 
uncertainty in best fit p is consistent with effects seen in 
1.3mm continuum maps of L1527 (Motte & Andre 2001). 

A standard Shu collapse model was not tested since the 
density distribution only flattens to p ~ 1.1 in a small 
region around the infall radius and is much steeper ev- 
erywhere else in the envelope. Zhou et al. (1996) were 
able to model molecular line emission with a Terebey-Shu- 
Cassen model (TSC, Terebey et al. 1984), a perturbation 
of the Shu collapse model to account for rotation. The 
one-dimensional averaged TSC proflle is similar to r~^'^ 
in the region inside the infall radius (5400 AU). A more 
detailed understanding of the density structure requires 
higher dimensional radiative transfer models to effectively 
model the asymmetric, flattened structure due to rotation 
and outflow cavities. 

4.6. L483 

L483 (IRAS 18148-0440) is another example of a Class 

protostar with large asymmetries observed in the submil- 
limeter continuum maps. The 850 and 450 fiia emission 
is extended in the northeast to southwest direction (Paper 
I), nearly perpendicular to the observed outflow direction 
(Parker 1988). Water maser emission has been detected 
towards L483 (Xiang & Turner 1992) as well as an associ- 
ated near- infrared nebula (Hodapp 1994). Molecular line 
studies identify this source as a possible collapse candi- 
date (Mardones et al. 1997) but with several confusing 
signatures due to infall and outflow (Park et al. 2001). 

The dust continuum emission from L483 is best fitted 
by a shallow {p = 1.2) power law with sisrf = 0.3 (Figure 
12). The model SED and 850 fiin radial profile match very 
well but the 450 /xm model profile is not steep enough in 
the outer half of the envelope. This may be partially due 
to uncertainties in the beam profile at 450 fiia, a variation 
in the dust opacity in the outer envelope, or a variation 
in the temperature structure due to geometrical effects 
compared to a one dimensional model. Once again, 0H5 
dust matches the observed SED well; however, the 160 /im 
point from Ladd et al. (1991) is much higher than the sur- 
rounding 100 fxm and 190 /xm points, perhaps indicating a 
larger calibration error than the published value. The in- 



ternal luminosity of the model is 13.0 Lq. This is the most 
luminous soiu'cx; in this sample of Class objects, mak- 
ing the contribution to the bolometric luminosity from the 
ISRF negligible compared to the internal luminosity, but 
the effects on the temperature profile are still substantial. 
The resulting temperature profile fiattens near the outer 
regions probed by SCUBA observations making a single 
temperature power law a poor approximation. Our model 
is consistent with 1.3 mm continuum models of Motte & 
Andre (2001), which assume a power law temperature dis- 
tribution, Td{r) oc r-«, {p = 1.2 ± 0.6 for g = 0.2 =F 0.2). 
The best fit model becomes optically thick around 30 /xm 
due to the flatter power law and higher luminosity internal 
source. 

Like L1527, the L483 profiles were not well fitted by a 
single power law in Paper I, but a power law density distri- 
bution does fit well when the beam, chopping, and Td(r) 
are correctly modeled. The best fit power law is close to 
p = 1, for both these sources with elongated dust emis- 
sion contours. Both L1527 and L483 have near-infrared 
nebulae associated with the core. However, both cores are 
not entirely similar since the extensions seen in the sub- 
millimeter continuum contours are in opposite directions 
with respect to the outflow axes. Geometrical projection 
effects are clearly important for these two sources. Tafalla 
at al. (2000) and Pezzuto et al. (2001) suggest this source 
is in the transitional phase between Class and Class I 
sources based on the observed properties of its outflow and 
far-infrared colors. While still a Class 0, it does have the 
largest Tbd (52K) of the sources we model here. While 
the envelope structure is clearly different from the more 
symmetrical cores B335 and B228, it is not clear whether 
the lower p indicates a more evolved state or an intially 
less spherical envelope. Higher dimensional dust modeling 
is required to fully understand the density structure of the 
outer envelope. 

No Shu collapse models were tested since the power law 
model was flatter than the Sim c;ollapse model. However, 
it is possible that a higher dimensional model that includes 
effects of rotation (TSC) may flt the intensity profile. 

5. DISCUSSION 
5.1. Density Distributions 

The properties of the best fit power law for all of the 

Class sources arc listed in Table 4. Power law models 
successfully fit the outer (r > 1000 AU) envelope struc- 
tures of Class protostars. Five of the Class sources 
are well fitted by steep power laws {p = 1.7 — 1.9, B335, 
B228, L723, IRAS03282-I-3035, L1448C), while the two 
cores with elongated emission contours are fitted by shal- 
low power laws {p = 1.1 — 1.2, L1527, L483). We quantify 
the degree of elongation in Table 4 by giving the aspect 
ratio (ratio of long axis to short axis) of the 20% contour. 
L1527 and L483 have the largest aspect ratios. Both L1527 
and L483 also have near-infrared nebulae, suggestive of as- 
pherical density distributions. The fiducial densities for all 
of the Class cores are consistent with the central densi- 
ties (uf ^ ric ^ 10^ cm~^) of the denser prc-protostellar 
cores modeled as Bonnor-Ebert spheres in Paper II. 

The best fit Shu77 models for B335, B228, L723, and 
IRAS03282-I-3035 have small infall radii (within the cen- 
tral SCUBA beam). The Shu77 models look very similar 
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to the best fit power laws, but generally do not fit as well. 
There is no strong evidence in the radial profiles for a break 
in slope indicative of an infall radius in the outer portion of 
the envelope probed by SCUBA. This result directly con- 
tradicts the molecular line modeling results towards B335 
(Choi et al. 1995), which strongly favors larger infall radii. 
However, the molecular line models do not take into ac- 
count chemical effects and abundance gradients, which can 
affect the shape of the blue asymmetry profile (Rawlings & 
Yates 2001). Unfortunately, the best fit power laws cannot 
be directly used in Monte Carlo molecular line modeling 
without an assumption about the velocity field along the 
line of sight. 

5.2. Mass Determinations 

The mass in the envelope within a 120" aperture was 
calculated using 
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The average mass is 2.6 ± 1.3 Mq for the seven Class 

best fit power laws, 3.5 times higher than the average pre- 
protostellar core mass (0.8±0.1 Mq) within the same aper- 
ture (cf. Paper II). The average mass within a fixed outer 
radius of 30,000 AU is 9.5 ± 5.1 Mq for Class sources 
compared to 2.7 ± 0.3 Mq for pre-protostellar cores. If 
these Class sources arc the next stage of evolution, they 
clearly have most of their mass in the envelope, as origi- 
nally proposed by Andre et al. (1993). 

In the absence of a realistic model for Td{r), masses are 
usually determined from submillimeter observations using 
an isothermal approximation. To facilitate mass determi- 
nations for sources without models, we explored the most 
suitable temperature to use. We calculated the isothermal 
dust temperature that yields the same mass as the detailed 
model based on the measured 850 /xm flux in a 120" beam, 
Sg^, using 
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The results are listed in Table 4. Most are 11-13 K, but 
the most luminous source, L483, has Tiso = 18 K. L1448C 
is excluded due to confusion from the northern sources 
in the 120" aperture. The mean value is 13.8 ± 2.4 K. 
The sources with p ~ 1 (L1527, L483) have higher T^so 
(15. OK and 18. OK respectively). A similar calculation can 
be made for the pre-protostellar cores modeled in Paper II 
using the mass of a Bonnor-Ebert sphere within the 120" 
aperture. For the three pre-protostellar cores, the average 
isothermal temperature is 11.1 ± 1.2 (12.4 K L1544, 10.9 
K L1689B, lO.OK L1512), slightly lower than for Class 
sources with an internal luminosity source. 

The model dust mass can be compared to the virial mass 
calculated using optically thin linewidths (e.g., N2H+ and 
H^^CO+; see Table 8 in Paper I). The virial mass esti- 
mates in Paper I assumed a constant density distribution 
in the envelope; however, the virial mass can be corrected 
for the best fit power law density distribution (cf. Bertoldi 
& McKee 1992) using 
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wlicrc 9ap is the FWHM size of the aperture in which the 
dust mass was determined, D is the distance, and {v^) is 
the 3-D velocity dispersion: 
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where Av is the FWHM linewidth of an optically thin line, 
/i is the mean molecular mass, and niamu is the molecular 
mass of the species whose linewidth is used. 

The mean of the ratio of the virial mass to the model 
dust mass in a 120" aperture is 2.1 ± 0.6. Uncertainties 
of up to a factor of 10 in the opacity exist between differ- 
ent dust models (Ossenkopf & Henning 1994). We regard 
this agreement as encouraging evidence that 0H5 opacities 
describe the dust in the cores well and that the virial theo- 
rem, properly applied, gives good mass estimates. This ra- 
tio is also consistent with the ratio of virial mass to model 
dust mass for a sample of deeply embedded high mass 
star forming cores associated with water masers (2.4 ± 1.4, 
Evans et al. 2002). A factor of 2 decrease in k;i^(850) from 
the 0H5 value would bring the dust model mass and virial 
mass into agreement on average. 

5.3. Spectral Indices 
The spectral index, (y.\fQ/gc,f), was calculated from the 
model fluxes at 450 and 850 /um in a 120" aperture 
(L1448C excluded due to confusion) using 

}^n-( Qvriod ( Qmod\ 
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log(850/450) " ^^^^ 

The model spectral index calculated within a 120" aper- 
ture agrees extremely well with the observed spectral index 

from Paper I; (^allQ/g^Q{model) /allQ/^^Q{observed)'j = 

0.96 ± 0.18. The model spectral indices are within the 
observed error bars; however, the total uncertainty on the 
450 /Ltm flux (^ 50%) makes the observed spectral index 
fairly uncertain. 0H5 opacities are also successful in re- 
producing the observed submillimeter spectral indices for 
Class sources. 

We also considered possible changes in the spectral in- 
dex as a function of impact parameter. To avoid cali- 
bration uncertainties we compare the normalized specific 
intensity, IlfS''^{b)/I^fS"'{b), for the best fit models to the 
observed ratios (Figure 13). The variations in the model 
ratio match those in the observed ratio. Apparent varia- 
tions in the spectral index can be removed by taking into 
account beam effects and a realistic temperature distribu- 
tion. For example, the first sidelobe at 450/xm results in 
a large increase in the ratio. The lower signal-to-noisc in 
outer annuli would mask any subtle variations in the dust 
opacity. L483 is a clear exception as the model specific in- 
tensity ratio does not decrease as fast as the observed ratio. 
The 450 /zm model is too flat at large radii. However, this 
discrepancy may be caused by our use of a spherical model 
on an aspherical source rather than an actual variation in 
the dust properties. There is thus little evidence for varia- 
tions of the dust opacity with radius. This result contrasts 
with some earlier work (e.g., Visser et al. 1998, Johnstone 
& Bally 1999) in regions forming more massive stars, but 
agrees with the conclusions of Chandler & Richer (2000), 
who were studying regions similar to those in this study. 
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It is extremely important to use a realistic beam profile 
and Trf(r) for modeling both the 850 and 450 /xm SCUBA 
maps. 

5.4. The ISRF and L,ni 

The contribution from the ISRF can be important to 
the heating of the region of the outer envelope probed by 
SCUBA observations of low luminosity cores. Scaling the 
ultraviolet to far-infrared portion of the ISRF by a fac- 
tor of 3 in either direction resulted in changes of the best 
fit power by factors of Ap = ±0.2. Every source profile 
except L1448C was better fitted by a lower strength of 
the ISRF, SisRF = 0.3, in agreement with models of pre- 
protostellar cores (Paper II). Modest shielding of the ISRF 
by surrounding gas in the molecular cloud surrounding the 
dense core could account for this. L1448C is forming in 
a much more crowded region and may be subjected to a 
stronger ISRF. 

Our models provide an estimate of the internal lumi- 
nosity of the central protostar. For all of our sources, the 
internal luminosity accounts for nearly all of the observed 
luminosity {{Lim/Lfj) = 0.95 ±0.12). While the heating 
from the ISRF is important for the temperature structure 
of the outer envelope, it does not contribute significantly 
to the overall observed luminosity of the Class sources 
we modeled. There is a wide range of internal luminosities 
modeled, from 1.0 Lq to 13.0 Lq, with an average inter- 
nal luminosity of 4.0 Lq. For each of the best fit models, 
the flux at 850 //m and the observed bolometric luminos- 
ity were fitted simultaneously. For sources with SEDs not 
observed shortward of 60 ^m, the model internal luminos- 
ity is a lower limit. All of the luminosities were calculated 
using the distances from Paper I. Since the distance to 
many globules arc uncertain to 50%, a true determination 
of the internal luminosity is uncertain to 100%; however, 
the model and observed internal luminosities are consis- 
tent for the distance adopted. 

Interestingly, the most luminous source, L483, has about 
the same envelope mass as the others. The greater lumi- 
nosity should reflect either a higher stellar mass or a higher 
accretion rate. If the former, it would imply a higher start- 
ing mass for the condensation; if the latter, it might be 
reflected in the linewidths. In fact, the linewidth of the 
N2II"'" line is similar to that of the other sources. Alter- 
natively, the accretion might be in a transient high state, 
similar to an FU Orionis event. 

5.5. Caveats and Future Work 

While our models take into account heating from a cen- 
tral source, heating from the ISRF, realistic beam effects, 
and simulated chopping, they cannot effectively model 
asymmetries seen in the dust continuum maps, flattening 
due to rotation or effects of magnetic fields, and clearing 
of material in outflow cavities. Five of the cores mod- 
eled (B335, B228, L723, IRAS03282+3035, L1448C) ap- 
pear sufficiently symmetric that the most important effect 
not included would be lower densities in outflow cavities. 
Harvey et al. (2001) found a strong asymmetry in the (H 
- K) colors along the outflow directions in B335, indica- 
tive of clearing of material along the outflow axes. The 
resulting best fit models included an outfiow opening an- 
gle of 35° to 45°. While no such asymmetry is observed 



in the submillimcter continuum maps of B335 (Paper I), 
extensions along the outflow directions (e.g., L1527, L723) 
and perpendicular to the outflow direction (e.g., L483) are 
observed in submillimcter maps of Class sources. The 
overall impact outflows have on the density structure of 
the envelope will become clearer with finer resolution and 
higher sensitivity. Future multi-dimensional modeling of 
the dust continuum emission should attempt to account 
for the effects of the outflow. 

We have neglected the emission from a disk. Chandler & 
Richer (2002) showed that the flux from a disk is negligible 
compared to the the total flux from the envelope of Class 
sources; however. Class sources may have a substantial 
fraction of emission from a compact component within the 
central beam. The interpretation of the density structure 
of the outer envelope may be strongly affected when a cen- 
trally normalized radial profile is used. Only a few Class 
sources have been observed at submillimcter wavelengths 
with interferometers (cf. Brown et al. 2000) resulting 
in few constraints on submillimcter disk fluxes. Observa- 
tions of 2.7mm continuum towards L1527 with OVRO and 
BIMA find a flux of ~ 40 mJy from a compact component 
(Terebcy et al. 1994, Shirley et al. unpublished obser- 
vation). Using a model for an active disk (T{r) ^ j.-0-5. 
Butner et al. 1994), we find an upper limit of 0.7 Jy at 850 
lj,m. In this scenario, the disk accounts for up to 70% of the 
flux within the central beam, decreasing p by Ap = —0.6. 
It is likely that the 7" aperture included flux from the en- 
velope, thereby overestimating the flux from a disk. If the 
disk has a steeper temperature power law (T(r) ^ r~° ''^), 
then the flux contribution drops to 0.1 Jy at 850 /xm and 
the change in the best fit power law model is Ap ~ —0.1. 
As another example, the Choi et al. (1996) Shu77 model 
fits the B335 radial profile when the disk fiux equals the 
envelope flux within the central beam (0.4 Jy at 850 /xm). 
However, BIMA millimeter observations of B335 do not 
support such a high disk flux (Shirley et al., unpublished 
observation) . Constraints on the disk flux and modeling of 
BIMA observations of the dust continuum towards Class 
protostars will be presented in a future paper. 

6. CONCLUSIONS 

We have modeled seven Class sources using single 
power law and Shu77 density distributions. Power law 
models suitably fit the 850 /im profiles and SED of all Class 
sources. Five sources with circular contours are best fit- 
ted by a steep value of the power law index p — 1.7 — 1.9 
(B335, B228, L723, IRAS03282+3035, L1448C), while two 
sources with aspherical emission are fitted by flatter power 
law indices p = 1.1 — 1.2 (L1527, L483). Uncertainties in 
the strength of the ISRF, sjsrf, and beam shape limit 
the accuracy in the power law index to ±0.2. 

The Shu77 model from Choi et al. (1995) does not fit the 
B335 radial profiles. Smaller infall radii are able to fit the 
profiles (B335, B228, L723, IRAS03282+3035, L1448C), 
but the infall radius is within the central 450 /um beam, 
effectively making the density distribution appear like a 
single power law throughout the region of the envelope 
probed by SCUBA. 

The average mass within a 120" aperture is 2.6M0 and 
is reasonably consistent with virial mass estimates and ob- 
served mass estimates from Paper I and models of the ini- 



tial conditions from Paper I. Wc find little evidence for 
variations in the dust opacity with radius. 0H5 dust re- 
produces the observed spectral index on average and pro- 
vides a good fit to many SEDs. In addition, it leads to 
masses determined from dust emission that arc consistent 
with virial masses to a factor of 2. Heating from the ISRF 
is very important for correctly interpreting the tempera- 
ture profile of the outer envelope of low mass star form- 
ing cores but does not significantly contribute to the total 
bolometric Imninosity. The dust models constrain the in- 
ternal luminosity of Class protostars, but distances to 
isolated cores remain the largest uncertainty in determin- 
ing accurate masses and luminosities. 

Outflow cavities and asymmetrical density distributions 
should be modeled using higher dimensional dust model- 
ing. In particular, L483 and L1527 have large asymmetries 
in the dust continuum emission that cannot be modeled 
with a one dimensional code. 

The presence of a disk within the central beam may af- 
fect the interpretation of the best fit density distribution in 
the outer envelope, decreasing p by as much as Ap = —0.6. 

We can now use a more realistic n(r) and Td(r) in Monte 
Carlo molecular line radiative transfer models to test infall 
models (which provide v{r)) and the origin of line asym- 
metries, to test predictions of chemical models, to inves- 
tigate amounts of depletion, and to improve estimates of 
the ionization fraction. 

Our primary conclusion is that a simple power law for 
the density distribution fits all of the Class sources that 
we have considered, while Shu77 models with substantial 
Tinf do not fit. Firmer conclusions await stronger con- 
straints on the submillimeter flux from a possible disk and 
modeling of interferometric observations of Class proto- 
stars. 
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Fig. 1. — Beam profiles used for modeling Class sources. The beam profile from January 1998, April 1998, and gaussian beams with the 
FWHMs reported in Paper I axe shown. Since CRL618 is a weak calibrator, the January 1998 beam was extrapolated beyond 34" using the 
April 1998 beam. 
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Fig. 2. — B335 Model SEDs for power laws with p=2.0 for various ro, ri, and Ku- Only the model SEDs are plotted since effects to the 
radial profiles are shown to be negligible. 
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Fig. 3. — B335 power Law models with p = 1.5, 2.0, and 2.5. The p = 2.0 model is the best fit to the SED and radial profiles. The model 
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Effects of changing the strength of the ISRF. Three p = 2 power law models of B335 axe shown with sjsrf = 0.3, 1.0, and 3.0. 
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Fig. 5. — The best fit power law model for B335. The parameters of the best fit were: p = 1.8, n/ = 1.5 x 10® cm Lint = 3.3L0, and 
SISRF = 0.3. 



18 



30 



o 



20 — 



10 — 







Tvl \l 1 1 1 

'a \ 1 
■•\ \ 

'•\ \ 
— \ 


1 M M 1 M M 
- 


'a \ 

"a \ 

"a \ 

'a \ 
"a \ 

"a 
"a 






\ \ 
■^ \^ ■ 

V . _ 
N 

^ 


Hnf 


6200 AU 




1000 AU 


Hnf 


1000 AU 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 



2.5 3 3.5 4 4.5 
Log r (AU) 

U I I I I I I I I I I I I I I I I I U 



— 



1 — 




850 /jLm 



I I I I I I I I I I I I I I I I I I I 



2.5 3 3.5 4 4.5 
Log b (AU) 



Q 



O 



— 



■1 



4 1111111111111 


1 1 1 1 


1 1 1 j- 


— SED 

- 






— 

: _M 






— .-^ > 






HI 






S|SRF 


1.0 ' 




S|SRF 


1.0 




~ S|SRF 


0.3 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 


1 1 1 r 



1 1.5 2 2.5 3 3.5 
Log X (microns) 







1 — 



J 1 M M M 1 M M 1 


MIL 














450 )Um 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


MM 



2.5 3 3.5 4 4.5 
Log b (AU) 
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Fig. 7. — The best fit power law model for B228. The parameters of the best fit were: p = 1.9, n/ = 1.2 x 10® cm , Li„t = l.OL©, and 
SISRF = 0.3. 
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Fig. 8. — The best fit power law model for L723. The parameters of the best fit were: p = 1.8, n/ = 1.8 x 10® cm , Li„t = 2.6L0, and 
SISRF = 0.3. 
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Fig. 9.— The best fit power law model for IRAS03282+3035. The parameters of the best fit were: p = 1.9, n/ = 1.9 x 10^ cm-3, 
Lint = 1.OL0, and sjsijF = 0.3. 
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Fig. 10. — The best fit power law model for L1448C. The parameters of the best fit were: p = 1.7, n/ = 2.6 x 10^ cm , Lint = 5.9L0, 
and sisRF = 1.0. 
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Fig. 11. — The best fit power law model for L1527. The parameters of the best fit were: p = 1.1, nj = 6.5 X 10^ Lint = 2.2L0, and 

SISRF = 0.3. 
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Fig. 12. — The best fit power law model for L483. The parameters of the best fit were: p = 1.2, n/ = 6.0 x 10® cm , Li„t = IS.OLq, and 
SISRF = 0.3. 
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Fig. 13. — The ratio of specific intensity, I^^q'^/Is^q'^ for the best fit power law models of isolated Class sources (solid line). Points are 
spaced at the half beam resolution at 850 /im. 
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Table 1 
Power Law Models 



Source 




Uf 




' o 






»S r c D IT 


Yorn 

A850 


A450 


^SED 


Notes 








fAU) 


fAU) 
















B335 


























2.0 


1.7 X 10^ 


60 


60 000 


0H5 


2.5 


1.0 


5.9 


3.9 


2.2 


Test Model 




2.0 


1.7 X 10*^ 


30 


60,000 


0H5 


2.5 


1.0 


6.8 


5.0 


2.1 






2.0 


1.7 X 10^ 


120 


60 000 


0H5 


2.5 


1.0 


5.4 


3.4 


2.7 






2.0 


1.2 X 10^ 


240 


60 000 


0H5 


2.7 


1.0 


4.4 


2.1 


5.9 






2.0 


1.2 X 10^ 


60 


120 000 


0H5 


2.5 


1.0 


5.0 


3.5 


5.7 






9 n 


9 n V 1 


fin 

LIU 


30,000 






1 n 


fi 


O.O 


u .o 






2.0 


3.0 X 10'^ 


60 


15,000 


0H5 


2.5 


1.0 


16.9 


10.8 


7.3 






1.8 


2.5 X lO*" 


60 


15,000 


0H5 


2.5 


1.0 


5.4 


3.9 


5.8 






1.5 


5.0 X 10^ 


60 


60 000 

\J\J ^\J\J\J 


0H5 


2.8 


1.0 


79.4 


34.5 


22.6 






2.5 


2.8 X 10^ 


60 


60 000 


0H5 


2.5 


1.0 


47.1 


19.7 


8.5 






2.0 


1.2 X 10^ 


60 


60 000 


0H5 


1.5 


3.0 


14.7 


9.5 


1.7 






2.2 


1.7 X 10^ 


60 


60 noo 


0H5 


1.5 


3.0 


5.6 


3.4 


5.6 






2.0 


2.3 X 10^ 


60 


60 000 


0H5 


3.0 


0.3 


12.2 


8.9 


4.2 






1.8 


1.5 X 10^ 


60 


60 000 


0H5 


3.3 


0.3 


3.8 


2.9 


19.6 


ijesu 




2.0 


6.0 X 10^ 


60 


60 000 

\J\J ^\J\J\J 


0H2 


4.4 


1.0 


6.2 


3.8 


101 




B228 


























2.1 


1.2 X 10^ 


60 


30 000 


0H5 


1.0 


1.0 


11.4 


8.6 


30.7 






1 Q 

-L • t7 


1.2 X 10^ 


fiO 


30 000 




1 

J- .u 


3 


8 8 


fi 4 


4 8 


Ktot 


L723 
























2.0 


2.0 X 10*^ 


60 


60,000 


0H5 


2.0 


1.0 


1.1 


1.3 


9.6 






1.8 


1.8 X 10*^ 


60 


60,000 


0H5 


2.6 


0.3 


0.8 


0.5 


3.4 


Best xlt 




1.8 


2.4 X 10^ 


60 


30,000 


0H5 


2.6 


0.3 


0.8 


0.6 


3.7 


IRAS03282'^ 
























2.1 


1.9 X 10^ 


60 


45,000 


0H5 


1.0 


1.0 


13.5 


3.5 


4.6 






1.9 


1.9 X lO'' 


60 


45,000 


0H5 


1.0 


0.3 


10.2 


1.2 


5.2 


Best xl 


L1448C 
























1.7 


2.6 X lO'' 


60 


45,000 


0H5 


5.9 


1.0 


4.7 


10.4 


2.0 


Best Xtot 




1.6 


2.8 X 10*^ 


60 


45,000 


0H5 


5.9 


0.3 


8.3 


9.4 


9.3 


L1527 


























1.5 


1.15 X 10^ 


60 


30,000 


0H5 


1.8 


1.0 


76.7 


39.6 


13.7 






1.0 


3.5 X 10^ 


60 


30,000 


0H5 


2.1 


1.0 


52.3 


60.3 


25.6 






1.1 


6.5 X 10^ 


60 


30,000 


0H5 


2.2 


0.3 


9.2 


5.4 


28.3 


Best Xr 


L483 


























1.3 


6.0 X lO'"^ 


60 


45,000 


0H5 


13.0 


1.0 


8.9 


5.1 


1.2 






1.2 


6.0 X 10^ 


60 


45,000 


0H5 


13.0 


0.3 


6.1 


5.2 


1.0 


Best Xtot 



"'XsED calculated using all flux points listed in Tables 5, 6, and 7 of Shirley et al. 2000 with A > 60/im. 
''IRAS03282+3035. 
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Table 2 

Sensitimty of Model Parameters" 



Variable Range Ap ALj„t 

OLq)_ 

ri 30 240 AU < 0.1 

To 30, 000 - 120, 000 AU < 0.1 0.2 

n/ 5 X 10^ - 3 X 10^ < 0.1 0.2 

Lint 1.5- 3.0 Lq <0.1 

0H5 vs. 0H2 < 0.1 1.9 

sisRF 0.3-3.0 ±0.2 1.5 

Pn{Of Jan. vs. Apr. ±0.1 

p 1.5-2.5 ... 0.3 



"Tested on B335. 

''The range of the variable tested. 

'^The beam shape. 



Table 3 
Shu77 Models 



Source 


'^infall 

(AU) 


tteff 

(km/s) 


V ■ 

(cm-3) 


ro 

(AU) 


(AU) 


Lint 


(Lg) 


X850 


X450 


XSED 


Notes 


B335 


























6200 


0.23 


60 


60,000 


0H5 


6.5 


1.0 


79.0 


53.7 


108 " 






6200 


0.26 


60 


60,000 


0H5 


5.0 


1.0 


82.1 


52.8 


62.2 






1000 


0.23 


60 


60,000 


0H5 


3.8 


1.0 


7.6 


3.6 


36.7 






500 


0.23 


60 


60,000 


0H5 


3.0 


1.0 


9.0 


3.5 


15.8 






1000 


0.23 


60 


60,000 


0H5 


4.5 


0.3 


2.3 


0.8 


66.8 ^ 


Best xl 




1000 


0.26 


60 


60,000 


0H5 


4.0 


0.3 


2.4 


1.0 


39.1 






6200 


0.23 


60 


60,000 


0H2 


6.5 


0.3 


39.6 


21.9 


59.0 ^ 






1000 


0.23 


60 


60,000 


0H2 


4.7 


0.3 


2.6 


1.1 


58.8 




B228 


























1000 


0.23 


60 


30,000 


0H5 


1.0 


0.3 


2.9 


15.6 


32.2 


Best Xtot 


L723 
























1000 


0.29 


60 


60,000 


0H5 


2.6 


0.3 


1.1 


0.6 


2.49 


Best xlt 


IRAS03282b 
























1000 


0.23 


60 


45,000 


0H5 


1.0 


0.3 


4.8 


13.3 


34.1 ^ 






1000 


0.26 


60 


45,000 


0H5 


1.0 


0.3 


3.9 


9.5 


18.7 


Best xlt 


L1448C 
























1500 


0.38 


60 


45,000 


0H5 


5.9 


1.0 


10.5 


3.0 


2.4 





"Unable to simultaneously match Ssso and L(,oj. 
''IRAS03282+3035. 
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Table 4 

Properties of Best Fit Models 



Source 


P 






Lobs 


Lint 


r mod 
^bol 


"450/850 


^1^0 


Aspecf^ 


J- ISO 


vir 


















Ratio 












'■'') 




(L,0 










(K) 




B335 


1.8 


1.5 X 


10*^ 


3.1(0.1) 


3.3 


3.1 


2.5 


2.6 


1.08 


13.2 


4.8 


B228 


1.9 


1.2 X 


10^ 


1.2(0.2) 


1.0 


1.1 


2.6 


0.8 


1.22 


12.5 


2.8 


L723 


1.8 


1.8 X 


10^ 


3.3(0.2) 


2.6 


3.4 


2.5 


3.9 


1.42 


12.4 


7.7 


IRAS03282'^ 


1.9 


1.9 X 


10^ 


1.2(0.3) 


1.0 


1.2 


2.4 


2.4 


1.03 


11.5 


3.6 


L1448C'= 


1.7 


2.6 X 


10^ 


6.0(0.5) 


5.9 


6.0 




4.5 








L1527 


1.1 


6.5 X 


10^ 


2.2(0.2) 


2.2 


2.2 


2.4 


1.6 


1.59 


15.0 


2.6 


L483 


1.2 


6.0 X 


10^ 


13(2) 


13.0 


12.9 


2.7 


2.3 


1.93 


18.0 


5.4 



"^Ratio of major to minor axis for the 20% peak contour. 
''IRAS03282+3035. 

^Quantities calculated using a 120" aperture not shown due to confusion from nearby sources. 



